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1. Introduction 

Lattice Monte Carlo calculations of the 
strong coupling constant a s = g 2 /4n have 
been the subject of considerable interest in 
the lattice community (see figure). Over the 
years different schemes have been devised 
and lattice calculations greatly improved, to 
the point that the recent world average is 
dominated by lattice results [1]. 

Values for a s are traditionally given at 
the Z-Boson mass in the MS scheme for Nf = 
5 fermion flavors. Alternatively, one may 
quote, e.g., the MS version of the dimen- 



world-averages w/o lattice (Bethke: 2000, 2012) 
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Figure 1: Lattice results for a a MS (Mz). Gray bars in the 
background mark the non-lattice averages from 2000 and 
2012 by Bethke (see, e.g., [1]) 
ter is often adopted when quoting Nf = 0,2 

lattice results. While schemes other than MS must be used on the lattice, a 1-loop calculation 



sionful QCD A parameter, A 



(Nf=5) 
MS 



The lat- 



suffices to yield A~ for different Nf in terms of the corresponding in these other schemes. 



MS 



(Nf) 

In 2007 [2] we introduced a method for determining using only Landau gauge gluon and 
ghost propagators. These 2-point functions are easily accessible on the lattice and a product of the 
gluon and ghost dressing functions, Zd and Zq, defines the strong coupling constant [3] 



< M (P) 



8l(c 



An 



Z D (p,a) Z%{p,a) (fora-r-0) 



(1.1) 



whose running at large scales (momenta) p is known up to 4-loop order [4]. In Eq. (1.1), #q is the 
lattice coupling at lattice spacing a and the suffix MM refers to the Minimal MOM scheme, the 
underlying renormalization scheme of this coupling [4]. Note that a^ M (p) is defined also beyond 
Landau gauge. The first four j3-function coefficients, j8o, . . . , ft, can be found for general covariant 
gauges in [4] . Meanwhile this coupling has also been used by a French-Spanish collaboration [5] 
to estimate A^ and a dimension-2 gluon condensate for different Nf. 

MS ° 1 



2. Simulation details 

Our data for a s MM is obtained on Nf = 0,2 and 2+1 SU(3) gauge field configurations. The 
quenched data is used to help quantify systematic errors, primarily due to finite volume and lattice 
discretization effects. Nf = configurations are thus generated (using the standard Wilson gauge 
action) for different lattice spacings between a/ro = 0.186 (j8 = 6.0) and a/r^ = 0.037 (jS = 7.2), 
and, simultaneously, for different physical volumes L 4 . We chose L/ro « 2.5,3.3 and 4.3, and the 
ro/a values from [6]. With our lattice spacings we cover a large range of momenta, reaching well 
into the perturbative regime (see below). This will allows us to fix '"oA^ free of assumptions on 
the source of deviations from pure 4-loop running at small momenta. 

The unquenched configurations are kindly provided by the QCDSF collaboration. For Nf = 2 
these are for the Wilson gauge action supplemented by clover-improved Wilson fermions, while 
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the Nf = 2+1 configurations are for a (tree-level) Symanzik improved gauge action and SLink 
fermions (see [7] for details). 

All gauge configurations are fixed to lattice Landau gauge using an iterative gauge-fixing al- 
gorithm. To guarantee high-precision the local violation of transversality is not allowed to exceed 
£ < 1(T 10 where, as usual, e = max* Ote Tr [(V M A^)(V /x A A - /i ) t ] and A xfl = 2^(U xfl - U^)\ bless . 

Standard techniques are also applied to calculate the gluon and ghost propagators. Their dress- 
ing functions are extracted, as usual, by using the exact tree-level structures of the respective lattice 
propagators. To reduce lattice discretization artifacts further, we consider data only for diagonal 
lattice momenta (strictest cylinder cut). Finally, data points for different j3 (and k) are brought onto 
the common scale r^p 2 , using the respective ro/a values. For Nf = we use the ro/a values from 
[6] and for Nf = 2 the (new) chirally-extrapolated values from QCDSF [8]. 1 For Nf = 2+1 we 
currently have data only for /3 = 5.50 for which ro/a rj 6.59 [9]. 

3. Quantifying the lattice corrections at large momenta 

Though the above standard tricks nicely smoothen the momentum dependence of our a s MM 
data, we cannot expect it to follow exactly the continuum 4-loop running, since, even away from 
lower momenta where nonperturbative effects (and further down finite-volume effects) become 
important, corrections due to the hypercubic lattice symmetry will be needed. For a high-precision 
analysis it is essential to quantify these corrections, which disappear only in the continuum limit. In 
particular for our Nf = 2 and 2+ 1 calculations a quantitative understanding of these corrections is 
important. For them the access to sufficiently small lattice spacings is still limited by the available 
configurations (a/ro rj 0.15 .. .0.12), and (small but systematic) deviations from the continuum 
4-loop running at large momenta can be clearly seen in the data (see Figs. 3 and 4 below). 

A complete removal of hypercubic lattice artifacts could be achieved using the so-called H(4) 
method, as in [5]. An alternative method, however, yields reduced statistical and systematic uncer- 
tainties. To motivate this approach, consider the gluon and ghost dressing functions, Zo and Z G , in 
(infinite-volume) lattice perturbation theory (henceforth "LPT"). At 1-loop order, e.g., these have 
the form 

Z D . G (ap) = l+g 2 Q [F D , G (a 2 p 2 ) +A D , G (ap)] °=° 1 + g 2 F D . G (a 2 p 2 ) (3.1) 
with Fo tG (a 2 p 2 ) = (c^j G log(a 2 p 2 ) +cfg G ) the part that survives in the continuum limit. The 1- 

D G D G 

loop coefficients c l0 ' and c n ' are known since the eighties [10]. For the hypercubic corrections, 
however, we are interested in Ad,g(^p)- These terms are non-zero for finite a, depend on the 
momentum ap and its direction (with ap^ = 2nk^/N tl and G (— N^/2,N^./2\] and represent 
the 1-loop contributions to the hypercubic lattice artifacts. For lattice hadron physics applications, 
e.g., it is common to calculate such correction terms and subtract them from the lattice data for some 
operators before determining the RGI scheme renormalization factors (see, e.g., [1 1]). Whether this 
removes sufficiently the bulk of lattice artifacts depends on the quantity being considered. 

For the gluon and ghost propagators, to the best of our knowledge, nothing is really known 
about Ad, g (cip). We therefore calculated them in LPT by evaluating (mostly by numerical integra- 
tion) the two Feynman diagrams (sunset and tadpole) for the ghost self-energy and another seven 

'Note that the new ro/a values have changed, this explains why ^qA^ is now larger than before [2]. 
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Figure 2: Hypercubic (1-loop) corrections, Aq (left) and A^ (right), versus a 2 p 2 . Aq is shown for on-axis 
(dash-dotted line) and diagonal momenta (solid line), while Ad only for diagonal momenta, but for Nf = 
(solid line) and Nf — 2 (dashed line) massless fermions. Note that at 1-loop, Aq is independent of Nf. 



for the gluon self-energy, all for the whole range of external momentum (ap). In particular the 
gluon self-energy diagrams containing 3- and 4-gluon vertices are quite cumbersome. 2 

The outcome of our 1-loop LPT calculation is shown in Fig. 2. The left panel shows Ao(ap) 
for two directions of ap, diagonal and on-axis. As expected, the hypercubic corrections grow with 
\ap\ and are smallest (but not at all zero) for diagonal ap. This is why the widely used cylinder cuts 
are so effective in removing a large fraction of the lattice artifacts from the data, though these are 
not removed completely. The right panel shows Ajj(ap) for Nf = and Nf = 2 massless fermions. 
For simplicity only diagonal momentum results are shown, the comparison of diagonal and on-axis 
momenta being similar to that seen for Ao(ap). 

Looking at Fig. 2 and the definition of a^ M one sees that, to leading order, the lattice cor- 
rections to a*"* (rap) are of the form a™ M (r p) + g 4 [A D (ap) + 2A G (ap)] + 0(g%) . Some of the 
corrections to the dressing functions will thus cancel in the a^ MM data at larger ap. 

After subtracting Ao(ap) and Ag(ap) corrections from our data for Zjj and Zq, we find, un- 
fortunately, only minor reductions in the size of lattice artifacts. Higher-loop corrections are thus 
needed, but calculations beyond 1-loop LPT are not feasible because of their complexity. Within 
the framework of Numerical Stochastic LPT this might be possible [14]. We have begun looking 
into this. 

In any case, from this exercise we know the remaining corrections should be of the leading 
order That is, a hypercubic Taylor expansion of the (1-loop) subtracted a s data for 

diagonal lattice momenta, denoted OC^c sub henceforth, should be of the form 

<C^{roP,ap) = CC™ M (r 2 p 2 )+g 6 [c 2 ■ (ap) 2 + c 4 ■ (ap) 4 + . . .] +O(g 8 ) (3.2) 

where C2 and C4 are constants. In fact, this leading-order expression describes the bulk of the lattice 
artifacts surprisingly well. We can fit our lattice data, CC^Hq mb (rop,ap), simultaneously for different 
/3 with the ansatz Eq. (3.2) with only three parameters: ci, C4 and toAmm (resp. roA^g), where the 
latter parametrizes the 4-loop running of a s MM (resp. oc^ s ) in the continuum. 

2 A.S. thanks Holger Perlt for his help and support of (many lines of) expressions for these vertices. 
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Figure 3: sub versus r oP 2 f° r = (top) and Nf = 2 (bottom). The gray band is a s MM {rip 2 ) running 
at 4-loop using r A^ = 0.637(32) and r A^ = 0.789(52) as input [12, 13]. The solid line (on this band) 
is aj MM (ro/? 2 ) as it comes out from a fit of 0!^ sub to Eq. (3.2) (no error band shown). For Nf = 2 all data 
shown has been used (simultaneously) for the fit. For Nf = we show data for only one physical volume 
for simplicity, but seven more sets entered the fit. The (colored) dashed curves, connecting the data points, 
refer to a^ M (rQp 2 ) plus the fitted lattice corrections (see Eq. (3.2)). Vertical lines mark the fit window: a 
dashed-dotted line marks the lower bound r^p 2 ^ = 800 (r^p 2 ^ = 600 for Nf = 2) and dotted lines the 
upper bound a 2 p 2 nax = 29. Open symbols are for points where a 2 p 2 > a 2 p 2 aax . 
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4. Results 

To illustrate how well ansatz (3.2) performs we show in Fig. 3 our (subtracted) a s MM data 
versus r\p 2 for Nf = (top) and Nf = 2 (bottom), together with the curves from a global fit. A solid 
(black) line refers to (tqP 2 ), the strong coupling in the continuum at 4-loop running, while 
the (colored) dashed lines refer to a, MM (rgp 2 ) + g% [c 2 • (ap) 2 + c 4 • (ap) 4 ] . Note that gl = 2N C /P 
varies with j3, but C2 and c\ are the same for all data sets of same Nf. For the fit shown in Fig. 3, 
the fitting window starts at r\p 2 = 800 (r^p 2 = 600 for Nf = 2) and ends at a 2 p 2 = 29. Other fit 
windows give similar results, but we see first deviations from 4-loop running below r\p 2 = 500. 
Due to the large lower bound finite- volume effects play no role. For the purpose of this illustration, 
we therefore fitted the data for all available volumes, and also did not distinguish the quark masses. 
This is suggested by the weak quark-mass (ff) dependence of our Nf = 2 data (the corresponding 
pion masses are m % 160. . .720MeV, see [8]). To not overload Fig. 3, the Nf = data is shown 
only for one physical volume. 

For a comparison, we also show the continuum coupling oc^ m (tqP 2 ) running at 4-loop (gray 
bands) using as input the recent values foA^ °' ) = 0.637(32) and H)A.j^ 2 ^ = 0.789(52) from the 
literature [12, 13]. Our fits fully agree with these values. 

As mentioned above, there is also a s 
data for Nf = 2 + 1. We are currently in 
the process of calculating a 4 MM on all the 
different sets provided by QCDSF to esti- 
mate t"oA~ . So far we have data only for 
one P (see Fig. 4) on rather course lattices 
(a/ro « 0. 15). But a first fit, similar to those 
discussed above, demonstrates already good 
agreement with the A^g from [15]. This is 
just a first fit attempt, and a more thorough 
study will follow. But our new approach 
of dealing with the lattice corrections (via 
Eq. (3.2) and restricting to data with diago- 
nal momenta) seems to work nicely. It will 

(Nf) 

allow us to fix >"oA^ , without the need to 
assume anything on the nature of the devi- 
ations from pure 4-loop running at smaller 
momenta (higher-loop corrections, dim-2 condensate, etc.) as was necessary in [5]. 

5. Summary 

We presented an update on our determination of rnA^ for different Nf. For this we use a 
renormalization scheme (Minimal MOM scheme in Landau gauge) that only requires lattice data 
for the gluon and ghost propagators in Landau gauge [4]. In this contribution we focused on a 
quantitative understanding of the hypercubic lattice artifacts at large momenta. This is important 
for a high-precision analysis, in particular for our unquenched calculations, for which the access to 
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Figure 4: Two sets of our Nf — 2 + 1 data and a first fit 
attempt similar to what was done for Fig. 3. 
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momenta above 10 GeV is limited by the available gauge configurations. To make maximal use of 
data at large lattice momenta (a 2 p 2 » 1) we determined the leading hypercubic lattice corrections 
to CC$£ mb . This was achieved by a calculation of the contribution at 1-loop order in LPT and 
secondly by a simultaneous fit of our lattice data (for diagonal momenta) to Eq. (3.2). It then very 
well describes the leading lattice corrections to sub with just three free parameters, one of 

(Nf) (Nt) 

which is roAj^ . Values for ^oA^ for Nf = and Nf = 2 (and perhaps even Nf = 2+1) will be 
given in a forthcoming article. There we will give also more details on our data and fits, and on the 
LPT calculations presented above. 

We are indebted to Holger Perlt for his help on getting the 1-loop LPT calculations right. This work 
is supported by the European Union under the Grant Agreements IRG 256594 and 249203. K.M. is sup- 
ported by NSERC (Canada) and L.v.S. by the Helmholtz International Center for FAIR. We thank the HLRN 
(Germany) for the generous support of computing time over the years. 
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